Frustrated spin-1/2 square lattice in the layered perovskite PbV03 
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We report on the magnetic properties of the layered perovskite PDVO3. The results of magnetic 
susceptibility and specific heat measurements as well as band structure calculations consistently 
suggest that the 5 = 1/2 square lattice of vanadium atoms in PbVC>3 is strongly frustrated due 
to next-nearest-neighbor antiferromagnetic interactions. The ratio of next-nearest-neighbor (J2) to 
nearest-neighbor (Ji) exchange integrals is estimated to be J2/J1 ~ 0.2 — 0.4. Thus, PbVOa is 
within or close to the critical region of the Ji — J2 frustrated square lattice. Supporting this, no 
sign of long-range magnetic ordering was found down to 1.8 K. 

PACS numbers: 75.47.Pq, 75.30.Et 
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The spin-1/2 square lattice provides a number of sim- 
ple and important models for theoretical physics. If one 
uses the Heisenberg Hamiltonian and considers nearest- 
neighbor (NN) interactions only, a well-known Heisen- 
berg square lattice (HSL) is formed. This model has 
been successfully applied to many transition-metal com- 
pounds (in particular, undoped high-T c superconductors) 
and much of its physics is now well understood^ A num- 
ber of new phenomena appear if next-nearest-neighbor 
(NNN) interactions are taken into account. If both NN 
(Ji) and NNN {J 2) interactions are antiferromagnetic, 
the spin lattice is frustrated since J\ and J2 tend to es- 
tablish different types of magnetic order. The properties 
of the system are controlled by the value of a = J2/J1 • If 
a is small, Neel order is favorable and the limit of HSL is 
realized. If a is large, columnar antiferromagnetic order 
is established. However, the main interest is attracted 
to the intermediate region that lies close to the quan- 
tum critical point at a c = 0.5. The nature of the ground 
state in this region is still disputed. Theoretical studies 
suggest different spin-liquid scenarios (e.g. RVB ground 
state) 

Most of the compounds realizing 5=1/2 square lat- 
tice are well described with HSL since NNN interactions 
are usually negligible. Recently two new systems with 
a ^> 1 (Li 2 VOX04, X = Si, Ge) were a subject of ex- 
tensive studies and revealed columnar antiferromagnetic 
order £J±il& However, no compounds in the critical region 
close to a c = 0.5 have been reported so far. 

Below we present magnetic properties of a novel com- 
pound, PbVO-3, 9 ' 10 that reveals 5 = 1/2 frustrated 
square lattice (FSL). Both experimental and computa- 
tional studies show that this compound lies close to the 
critical region of FSL and does not undergo long-range 
magnetic ordering down to 1.8 K. 

PbVC>3 adopts a layered perovskite-type structure 
(space group P4mm, a = 3.8001 A, c = 4.6703 A) 
shown in Fig. [T] This structure combines the absence 
of an inversion center with a magnetic V + cation, there- 
fore PbVC>3 attracts considerable attention as a possi- 
ble multiferroic i 11 i 12 ' 13 i 14 Magnetic properties of PbVC>3 





FIG. 1: (Color online) Perspective view (left panel) of the 
crystal structure of PbVOs and projection along the [001] 
direction (right panel). VO5 square pyramids share common 
corners, lead atoms are indicated by large spheres 



have not been reported so far. 

Polycrystalline samples of PbVC>3 were prepared by a 
high-temperature high-pressure technique in a belt-type 
apparatus. Stoichiometric mixtures of PbO and VO2 
were placed into gold capsules and treated at 950 °C for 2 
hours under a pressure of 5 GPa. The phase composition 
of the prepared samples was checked by XRD. 

The magnetic susceptibility of PbVC>3 was measured 
between 1.8 and 400 K in the fields [ioH of 0.1, 1, and 5 
T using a Quantum Design SQUID magnetometer. The 
specific heat was studied in the temperature range from 
1.8 to 270 K with a commercial PPMS. 

The study of the magnetic properties of PbVC>3 is com- 
plicated by the presence of magnetic impurities in the 
samples under investigation. All of the measured sus- 
ceptibility curves revealed anomalies at 90 and/or 350 K 
corresponding to PbVgOn and VO2 respectively^ al- 
though these impurities were not necessarily visible in 
XRD patterns. Fig. [5] presents one of the best x(T) 
curves (no anomaly at 90 K) below 300 K, since the high- 
temperature part is affected by VO2. The features of this 
curve are typical for all PbVC>3 samples within a variation 
of the temperature-independent background from sample 
to sample due to different amount of VO2. These features 
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FIG. 2: (Color online) Magnetic susceptibility of PbVOs mea- 
sured at 5 T, arrows denote the curves obtained at heating 
and cooling. The dashed and solid curves in the inset show 
the fits by high-temperature series expansion for HSL model 
and FSL model with J2/J1 = 0.38, respectively. 



are: i) a very broad maximum near 190 — 200 K; ii) the 
difference between the zfc and fc curves below about 50 K 
with small humps near 43 K (or heating/cooling hystere- 
sis visible at Fig. [2]). One should note that the anoma- 
lies at 45 — 50 K may appear due to the influence of 
trapped oxygen that undergoes condensation at about 
50 K. These anomalies are especially strong if the in- 
trinsic signal is very small (as in the case of our study). 
Nevertheless, the susceptibility data alone do not allow 
to decide unambiguously whether the anomaly at 50 K 
is intrinsic to PbV03 (e.g. indicating a phase transition) 
or not. Therefore, we turn to other experimental data in 
order to check the extrinsic nature of this anomaly. 

The specific heat curve (Fig. \5§ is smooth between 1.8 
and 270 K and suggests the absence of phase transitions 
in PbVC>3 in this temperature range. The conclusion 
is supported by thermal expansion 1 ^ and resistivity^*^ 
data. Finally, a neutron diffraction study at 1.5 K shows 
the absence of long-range magnetic ordering^ Thus, we 
conclude that the susceptibility anomaly at 50 K has ex- 
trinsic nature. 

Square-pyramidal coordination of vanadium is known 
to give rise to a non-degenerate d xy ground state for V +4 
(see Refs. IMT7I). The d xy orbitals lie parallel to the V-0 
layers, therefore one may expect that PbVC>3 is a 2D 
spin system. Indeed, the broad maximum in the sus- 
ceptibility curve (Fig. \2§ is typical for 2D spin systems 
while the upturn at low temperatures is usually ascribed 
to the paramagnetic contribution of impurities and de- 
fects. Unfortunately most of the regions of the curve 
are unsuitable for fitting since the low-temperature part 
is affected by impurity contributions while the metal- 
insulator transition of VO2 prevents us from using the 
data above 300 K.— Thus, only the narrow region be- 
tween the maximum (T max = 190 — 200 K) and 300 K is 
appropriate for the fitting. Since the leading exchange in- 
tegral J is comparable to T max , high-temperature series 




FIG. 3: (Color online) Specific heat of PbV0 3 . The inset 
shows the comparison with the diamagnetic reference com- 
pound PbTi0 3 . 



expansions (HTSE) are applicable in this region. 

We use two types of HTSE corresponding to HSL 
and FSL models, respectively^^ A temperature- 
independent contribution was added to the series as an 
adjustable parameter (\o) in order to account for the 
contribution of VO2 as well as other (diamagnetic or 
van Vleck) contributions of this type. Both models re- 
sulted in fits of similar quality with reasonable values of 
Xo ~ 2.5 • 10~ 4 emu/mol. We find J = 190 K for HSL 
and Jt = 203 K, a = J2/J1 = 0.38 for FSL. The g value 
was fixed at g = 2 in order to get a more stable fit. Ba- 
sically, we see that the susceptibility curve is consistent 
with both scenarios - HSL or FSL. 

Specific heat data may provide additional quantitative 
information about the spin system of PbVC>3. However, 
one has to estimate and subtract the phonon part of the 
specific heat in order to reveal the magnetic contribution. 
The leading exchange integral in PbVC>3 is about 200 K 
therefore the magnetic contribution does not fall to zero 
even at high temperatures and can not be separated from 
the phonon part by a simple Debye-fit. A reference dia- 
magnetic compound has to be found in order to give a 
reliable estimate of the phonon contribution. 

We tried to use PbTiC>3 as a non-magnetic reference. 
The crystal structures of PbVC>3 and PbTiC>3 look sim- 
ilar but vanadium and titanium adopt different coordi- 
nation (square pyramid and octahedron, respectively). 
This difference seems to be crucial: the c p (T) curves have 
crossings at low temperatures (see the inset of Fig. [3]) in- 
dicating quite different phonon spectra. We are not aware 
of any non-magnetic compound with layered perovskite- 
type structure perfectly matching that of PbV03. There- 
fore, presently we can not estimate the magnetic contri- 
bution to the specific heat of PbV03. 

Experimental data provide rather limited informa- 
tion about the magnetic properties of PbVC>3, there- 
fore we turn to computational techniques. Band struc- 
ture calculations are known as a useful tool for esti- 
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FIG. 4: (Color online) LDA band structure of PbVC>3 near 
the Fermi level. Dots show the contribution of V 2>d xy states, 
thick solid line presents the fit of the tight-binding model. 
The inset shows the density of states for PbVOa near the 
Fermi level. 



mating exchange integrals and understanding the prop- 
erties of low-dimensional spin systems. In particular, 
computational results were helpful in the study of FSL 
compounds Li2VOX04 and provided the first realistic 
estimate 5 a ^> 1 (supported by later neutron experi- 
ments, see Ref. 0) in contrast to the early experimental 
reports^ 

Scalar relativistic band-structure calculations were 
performed using the full-potential local-orbital scheme 
(FPLO, version 5.00-19)22, and the parametrization of 
Perdew and Wang for the exchange and correlation 
potential^ A k mesh of 1152 points in the Brillouin zone 
(224 in the irreducible part) was used. 

We start with the LDA band structure (Fig. 2]). The 
highest occupied band reveals a significant contribution 
of V d xy orbital consistent with the square-pyramidal 
coordination of V +4 (see Refs. IMT7TI. A simple one- 
band tight-binding model is fitted to this band. We find 
t x = 0.132 eV and t 2 = 0.077 eV for NN and NNN 
hoppings, respectively (see Fig. [TJ. Long-range in-layer 
hoppings do not exceed 0.004 eV and hence may be ne- 
glected. The inter-layer hopping is also extremely small 
(t± = 0.0012 eV) suggesting a strongly two-dimensional 
spin system in PbVC>3. The t values are used to estimate 
antiferromagnetic contributions to exchange integrals as 
jAfm _ 4^|^jj eff (U e s is an effective on-site Coulomb re- 
pulsion). Assuming U c s = 4 eV we find J\ = 203 K, 
J± w 0.01 K, a = t\ji\ — 0.34 in a perfect agreement 
with the FSL fit of the susceptibility data. Note that the 
a value is a direct result of the tight-binding fit and does 
not depend on {7 e ff. 

Thus, the tight-binding model supports the FSL sce- 
nario for PbVC>3. However, this model assumes that 
all the interactions in the system are antiferromagnetic. 
Sometimes it is not the case and one may calculate total 
energies for different spin states in order to check pos- 
sible ferromagnetic contributions to the exchange inte- 



TABLE I: The total energies of different spin states of PbVC>3 
(the values are given in meV relative to the energy of FM 
state). 
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grals. LSDA calculations for several simple spin states 
have been reported ) 11 ' 12 and we list these results in the 
first line of Table [J 

LSDA results for transition metal compounds are 
known to suffer from an unreasonable neglect of corre- 
lation effects. Below we show that the relative values of 
total energy are considerably changed if one takes into ac- 
count strong Coulomb correlation for vanadium 3c? shell. 
The most simple way to introduce such correlation within 
DFT is provided by LSDA+[7 technique^ We apply sev- 
eral physically reasonable values of U (further named Ud 
in order to distinguish it from J7 e ff) and fix J = 1 eV. 

Table [J reveals quantitative rather than qualitative 
dependence of the LSDA+J7 results on the Ud value. 
The energy differences (£a-Afm — -Efm) and (-E-c-Afm — 
-Eg-afm) corresponding to inter-layer coupling J± are 
now comparable and small in contrast to the LSDA re- 
sults. The sign of the inter-layer interaction still is not 
clear but the absolute value of J± has the order of units 
of K. The in-layer interaction exceeds Jj_ at least by an 
order of magnitude. Thus, the spin system of PbV03 
is two-dimensional consistent with the crystal structure 
and the tight-binding results. 

Now we use LSDA+J7 to estimate in-layer interactions 
J i and J2. The results listed in Table HT1 are in a rea- 
sonable agreement with the experimental data and the 
tight-binding model, although the a value is somewhat 
lower than 0.34 - 0.38. Nevertheless, LSDA+C7 indicates 
considerable NNN interaction. Thus, band structure cal- 
culations strongly support the FSL rather than the HSL 
scenario. 

Additional evidence for the frustration in PbV03 is 
found if one considers the presence of long-range mag- 
netic ordering in this compound. The HSL tends to long- 
range order even at very weak inter-layer coupling (for in- 
stance, Neel temperatures Tn of undoped high-T c super- 
conductors have the order of hundreds of K, see Ref. I23T) . 
If J± in PbV03 is as small as 0.01 K (tight-binding result) 
T N ~ 0.2Ji w 40 K^i The phase transition at 40 K can 
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hardly be missed while analyzing experimental results. 
However, frustration effectively reduces Tn or even pre- 
vents the system from long-range magnetic ordering at 
all. 

Summarizing, all the experimental data available for 
bulk samples of PbVC>3 do not provide evidence for long- 
range magnetic ordering down to 1.8 K. Thermodynamic 
data (magnetic susceptibility, specific heat, thermal ex- 
pansion, and resistivity) do not show any indications for 
intrinsic phase transitions, while low-temperature neu- 
tron diffraction reveals the absence of long-range spin 
ordering at 1.5 K. 

Note that the recently reported thin films of 
PbVO ^ 13 ^ 14 show a different magnetic behavior with a 
possible magnetic ordering at 100-140 K. However, bulk 
phase and thin film may have distinct properties due to 
slightly different structures and strain effects. Detailed 
structural information for thin films of PbV03 is not 
available, but the reported difference of the c parameter 
of the unit cell (4.67 and 5.02 A in bulk and thin film, re- 
spectively) suggests considerable change of the structure 
in thin film as compared to bulk solid^ 

Thus, we conclude that magnetic frustration of the 
square lattice is a crucial feature of PbVC>3. If a > a c i 
[a cl = 0.34 (Ref. [3) or 0.24 (Ref. S)] long-range ordering 
is suppressed and some kind of spin-liquid ground state 
is formed. Unfortunately, we can not give a more precise 
estimate of the a value for PbVC>3, hence we can not con- 
clude if PbVC>3 falls into the spin-liquid region or is just 
close to its boundary. Anyway, the frustration in PbVC>3 



is strong enough to suppress Tn considerably (at least by 
an order of magnitude as compared to the HSL estimate) 
or provide a disordered (spin-liquid or spin-glass) ground 
state. 

Summarizing, this study provides strong evidence that 
PbVC>3 is the first example of S = 1/2 square lattice sys- 
tem lying within or close to the critical region of the FSL 
phase diagram and lacking for long-range magnetic order- 
ing down to low temperatures (Tn/J < 0.01). A further 
study of the ground state and low-temperature proper- 
ties of this system is of high interest. Unfortunately, we 
have to point out considerable difficulties in the synthe- 
sis of bulk single-phase samples (and, moreover, single 
crystals) of PbVC>3. The recent study of epitaxial thin 
films of PbVO j 13 ' 14 provides an alternative route for the 
preparation of this interesting compound but thin films 
and bulk solids may have strikingly different properties 
as we have mentioned above. 

In conclusion, we have shown that PbVC>3 reveals sig- 
nificant frustration of the square lattice that prevents this 
compound from long-range spin ordering down to 1.8 K. 
PbVC>3 lies very close to the critical region of the FSL 
phase diagram, and it may provide the first realization 
of disordered ground state for S = 1/2 square lattice. 
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